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1.  INTRODUCTION 

Certain  consistencies  have  bqen  shown  to  exist  among  the  anomalies 
of  short  period  teleseismic  P-phase  amplitudes  when  those  anomalies  are 
grouped  according  to  geographical  regions,  in  this  report  we  will  present 
evidence  that  the  source  path  strongly  influences  the  observed  amplitude 
anomalies . 

2.  PROCEDURE 

Eight  earthquakes  from  the  same  geographical  area,  viz.,  the  Fiji 
Islands,  were  selected  for  this  investigation.  All  eight  events  occur¬ 
red  at  143°  Az  with  respect  to  the  center  of  the  LASA  and  ranged  from 
9,500  km.  to  10,500  km.  distance.  These  events  were  ordered  in  relation 
to  their  distance  from  the  LASA  and  were  partitioned  into  pairs.  The 
result  of  this  grouping  is  a  set  of  four  pairs  of  events  such  that  the 
elements  within  each  pair  have  epicenters  within  eight  kilometers  of  each 
other.  We  note  here  that  although  these  epicenters  are  close  together 
the  focuses  of  the  events  have  a  considerably  greater  range.  More  will 
be  said  later  regarding  this  aspect.  We  now  list  the  paired  events  and 
some  of  their  identifying  characteristics. 

TABLE  1 


Pair 

1 

Pair 

2 

Pair  3 

Pair  4 

Event  No. 

152 

253 

219 

254 

197 

247 

171 

342 

Azimuth 

243° 

243° 

243° 

243° 

243° 

243° 

243° 

243° 

Distance ,km. 

9569 

9565 

9536 

9584 

10170 

10162 

10423 

10422 

Depth ,  km . 

33 

35 

127 

33 

338 

320 

525 

511 

Magnitude 

5.0 

5.1 

4.9 

5.2 

4.8 

5.5 

5.1 

5.4 

1 


Each  pair  consists  of  two  independent  earthquakes  in  the 
sense  that  all  were  separated  in  ^ime  by  two  to  ten  weeks.  Each 
of  the  pairs  was  investigated  for  the  presence  of  correlation 
between  member  events  in  the  following  manner.  The  peak  to  peak 
P-phases  amplitudes  recorded  by  the  center  elements  of  the  sub¬ 
arrays  were  determined  for  each  event .  In  this  way,  an  ordered 
set  of  amplitudes  indexed  by  the  recording  instruments  was  as¬ 
sociated  with  each  event.  Thus,  in  order  to  detect  the  presence 
of  correlation  between  any  two  events  one  merely  tests  for 
correlation  between  the  two  sets  of  amplitudes  generated  by  the 
respective  events.  But  before  taking  this  step,  it  is  important 
to  note  that  an  accurate  interpretation  of  the  correlation  coef¬ 
ficient  is,  in  general,  possible  "only  when  the  underlying  popu- 

(2) 

lation  is  normally  distributed."  It  has  been  demonstrated 

that  a  loqnormal  distribution  describes  the  behavior  of  amplitude 

(1) 

anomalies  over  LASA  with  remarkable  precision.  Hence,  by 
using  the  logarithms  of  the  peak  to  peak  amplitudes,  we  have 
characterized  each  event  with  a  normally  distributed  set  of 
measurements  and  can  now  apply  the  correlation  theory.  The 
estimate  of  the  coefficient  of  correlation  is  defined  by 

c  « 

r  =  -  ■  Cov(x»Y)  — .  The  following  table  presents  the  calculated 
J  Var  X  •  Var  Y 

estimate  of  the  coefficient  of  correlation  along  with  the  5 %  and 
1%  points  for  the  equal-tails  test  of  the  hypothesis. 


TABLE  2 


Pair  1 

Pair  2 

Pair  3 

Pair  4 

r 

0.613 

0.257 

0.409 

0.624  j 

5 %  Point 

0.482 

0.154 

0.532 

0.602 

1%  Point 

0.606 

0.641 

0.661 

0.735 

Conclusion 

Correlated 

Correlated 

Correlated 

Incon¬ 

clusive 

One  pair  appears  to  be  highly  correlated,  two  to  be  uncor¬ 
related,  and  one  to  be  a  borderline  case.  Ranking  these  pairs 
from  highest  to  lowest  degree  of  correlation  yields  1,  4,  3.  2. 

If  we  also  rank  the  pairs  from  those  with  earthquakes  which  occur¬ 
red  closest  together  to  those  with  earthquakes  farthest  apart,  we 
get  — — (4.v5  — 4 — (14.1  km),  3  (19.7  km),  2  (94  km)  ;  that 

is,  the  pairs  have  the  same  ordering  in  both  cases.  The  events 
which  occurred  closest  together  had  the  highest  correlation  and 
those  vhich  occurred  farthest  apart  were  least  correlated.  This 
relationship  is  better  illustrated  in  the  accompanying  drawing. 
Figure  1A,  which  depicts  the  distance  between  the  focuses  of  the 
paired  events  and  corresponding  estimates,  r. 


To  test  further  the  apparent  relationship  between  degree  of 
correlation  and  physical  proximity  of  events,  all  remaining  twenty- 
four  possible  combinations  of  paired  events  were  tested  for  cor¬ 
relation.  Since  the  separation  between  events  in  these  combinat¬ 
ions  (except  for  four  pairs)  is  much  greater  th?n  that  of  the  four 
pairs  discussed  above,  one  would  expect  that  there  would  be  an 
absence  of  correlation  in  these  cases.  Table  3  in  Appendix  A 
lists  the  computed  correlation  coefficients. 


The  two  pairs  marked  by  an  asterisk  in  Table  3  are  consistent 
with  the  results  illustrated  in  Figure  1A  and  in  Table  2.  Speci¬ 
fically,  we  find  (152,254)  has  r  =  0.613  and  separation  of  15  km 
and  *(253,254)  has  r  =  0.521  and  separation  of  19.1  km.  There 
still  remain  five  pairs,  three  borderline  and  two  correlated, 
which  appear  inconsistent  with  Figure  1A  and  Table  2.  This  may 
suggest  that  these  events  are,  in  some  sense,  "in  phase''.  But 
also  we  note  that  there  is  one  chance  in  twenty  (at  the  5%  point) 
that  we  may  compute  a  coefficient  which  implies  correlation  when  none 
exists.  Hence  we  should  "expect"  that  at  least  one  of  the  three 
borderline  cas<»s  is  not  really  correlated. 


Table  4  (Appendix  A)  shows  the  approximate  separation  between 
all  of  the  event  combinations  listed  in  Table  3.  A  graph  of. cor¬ 
relation  coefficients  versus  separation  of  focutee  ie  presented  in 
Figure  IB  for  the  four  pairs  in  Table  1  and  the  four  pairs  in  Table 
4  whose  focus  separation  was  less  than  100  km. 


3 .  ADDITIONAL  EVENTS 

Similar  tests  were  performed  on  three  teleseismic  events  which 
occurred  near  the  surface.  It  is  believed  that  these  three  origina¬ 
ted  at  the  same  location.  All  possible  combinations  of  pairs  were 
tested  and  in  each  instance  a  high  positive  correlation  was  found 
to  exist.  It  is  fairly  certain  that  these  events  did  not  all  occur 
at  exactly  the  same  azimuth  and  distance,  and  this  may  account  for 
the  variation  in  the  computed  coefficients.  The  coefficients  were 
computed  to  be  0.900,  0.671,  and  0.837  for  the  three  combinations 
of  pairs.  Figures  2,  3,  and  4  are  plots  of  the  logarithms  of  the 
subarrays  center  elements  for  ell  combinations  of  the  three  eventB. 


4.  ERRORS  OF  MEASUREMENT 

Hald^3^  devotes  a  section  to  the  influence  of  errors  of  measure 
ment  on  the  correlation  coefficient.  Briefly,  assume  (x^  ,  x2)  are 
two  true  properties  whose  relationship  is  being  studied.  If  the 
errors  of  measurements  are  denoted  (v^  ,  v^)  the  observations  are 
then  (y  ,  y2)  =  (x1  +  v,  ,  x2  +  V2**  Assumin9  these  errors  are 
stochastically  independent  and  independent  oE  (x^  ,  x^) ,  the  coef¬ 
ficients  p(v^  ,  v2)  =  p(x1  ,  v^)  =  p(xx  ,  v2)  =  p(*2  '  vi^  =  p^x2  ' 

=  0.  Assuming  the  means  of  v^  ,  ,  x^  ,  x2  are  zero,  we  have 


2  2 

Var (y . )  =  Var(x.)  +  Varfv^)  CTX  +  av  <■ 
11  i  i 

of  Covariance, 

Cov(y1  ,  y2)  =  E  f(x1  +  vx)  •  (x2  +  v2)] 


i  =  1,2  and  by  definition 
-  E(x1  +  vx)  •  E(x2  +  v2) . 


-  4  - 


Since  E(x^)  =  E{v^)  =  0,  i  =  1,2,  we  have 

Cov(y1  ,  y2)  =  e[  (x1  +  v^  -  (x2  +  v2>  j 

=  e[xiX2  +  +  x^  +  Vlv2  j 

=  Etx^)  +  Etx^)  +  E{x2v1)  +  E(v  v  ) 

but  by  the  assumption  of  independence,  we  note  that 

Etx^)  =  Efx^  *  E{v2)  =  0-0  =  0; 

E(x2Vi>  =  E(x2)  -  E(vx)  =  0-0  =  0; 

E{ViV2)  =  E{v1)  *  e(v2)  =  °*°  “  °* 

and  this  leaves 

Cov(y^  ,  y2)  =  E{x1x2).  But  recalling  that 
Etx^  =  0  and  E(x2>  =  0,  we  get 
Cov(y1  ,  y2)  =  E{x1x2)  -  0-0 

=  Efx^)  -  E(x1)  E(x2) 

Cov (yx  ,  y2)  =  Cov(x1  ,  x2) 


Now  by  definition,  p  (x^  ,  x2)  = 


Cov(x1  ,  x2) 


or  (Cov^  ,  x2)  =  [P(xx  ,  x2)j  •  (a*  )  -  (ux  )  and  finally 

l  *2 


<’'Ov(y1  ,  y2)  =  [p{Xl  ,  x2)] 


The  definition  given  by  (I)  then  says 


'  .  x.  • :  v  ■'  ■■■rr 


:v :  ?..»;••• 

vv^^sgip 

.’■’  v  v«  r  ::_■  j»,«v-  «»' s;- -c 


This  says  the  correlation  coefficient  of  the  true  properties 
plus  errors  of  measurement  is  always  less  than  the  correlation 
coefficient  of  the*  properties  themselves.  Equation  (II)  shows 
that  relatively  large  errors  of  measurement  disguise  existing 
correlations.  Therefore,  we  may  generally  say  that  the  computed 
coefficient  is  always  going  to  be  smaller  than  true  p. 

5.  DENSITY  FUNCTION  OF  r 

Since  we  have  been  working  with  relatively  small  samples, 
we  were  interested  in  the  behavior  of  the  density  function  for  r 
when  n  is  small.  R.  A.  Fisher  is  credited  with  discovering  the 


distribution  of  the  correlation  coefficient. 


This  rho  distri¬ 


bution  is  a  function  of  n,  the  sample  size  and  p,  the  correlation 
coefficient.  One  expression  for  the  density  function  is 


- 


n-1  r 


o  (1-pry)" 


V  l-x2 


The  integral  term  can  be  simplified  by  setting  x  -  Sin  0,  dx  = 

Cos  OdD  and  correcting  the  limits  to  0,n/2  I.  The  integral  becomes 
n/2  n  2 

P  Sin  0  CobQ  dQ _ 

o  {1-p.TSin  a)  "-1 

.  .  tt/ 

Noting  the  identity  we  get  ^  2  0.„n-2  n 

J 


2  2 
1  -  Sin  0  =  Cos  0 


Sin11  2  0  d0 


-*■  W4.II  W  WO  O  «  /  .  n  .  II  X 

o  (1-prSmO) 

Not  finding  a  closed  form  solution  for  this  integral,  a  numerical 
solution  using  Romberg  integration  was  devised  by  Dr.  C.  S.  Duris. 

When  n  =  2,  it  is  obvious  from  (II)  that  f2(r)  =  C  and  there¬ 
fore  the  coefficient  of  correlation  is  necessarily  +  1.  Graphs  for 
n  =  c,4,6  and  50  are  presented  Figures  5,6,7  and  8  in  Appendix  C 
along  with  a  listing  of  the  computer  program. 
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TABLE  3 
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Paired  Events 


i  ft  >M  $1 
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Borderline  (15  km) 

<- 

Borderline  (19.1  km) 
No  Correlation 

H  ii 


W  0, 


Borderline 


IB 


5%  Point 

1%  Point 

0.482 

0.666 

0.532 

0.754 

0.497 

0.497 

0.532 

0.514 

0.482 

0.576 

0.497 

0.623 

0.497 

0.623 

0.553 

0.497 

0.482 

0.666 

0.514 

0.641 

0.497 

0.623 

0.482 

0.606 
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